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The present invention relates to a method of measuring the 
dilution of Phase modulated spins as well as to a magnetic reso- 
nance imaging device which can be used for performing thxs 

m6th0 The present invention finds particular application in regard 
to medical diagnostic magnetic resonance imaging, but it is to be. 
appreciated that the present invention also finds application m 
n^netic resonance spectroscopy and magnetic resonance imaging 

^ ""g"sues water protons are often compartmental- 
ized in terms of their molecular mobility. Commonly, the two-pool 
m0 del is used to explain macroscopic relaxation properties re- 
sulting from this compartmentalization: One pool, usually the 
largest in size, is associated to the -free" mobile bulk water, 
« a second pool is associated to motional -triced pro- 
tons bound to macromolecules such as proteins or ipids In Mag 
netic Resonance in Medicine- (Vol. 10, P. 135-44 (1989)) S. D. 
Wolff and R . S . Balaban first coined the terms free and re- 
stricted proton pools to describe these different pools. Most im- 
portantly both pools are coupled by an exchange of magnetization 
TaTemlcal exchange and dipolar coupling. This phenomenon, com- 
monly termed magnetization transfer (MT) , is described in U S. 
Pat. No. 5,050,609 and how to make use of it in magnetic reso- 
nance imaging is reviewed in "Magnetic Resonance Quarterly" (Vol. 

8 P 116-37 (1992)). . • , 

in medical diagnostic imaging the quantity of the restricted 
< ma cromolecular) proton pool is highly desirable as it is ex- 
pected that this pool directly represents tissue structure^ 
tissue integrity. However, depicting of the restricted proton ^ 
pool is not possible because in most tissues the transverse mag- 
netization of these protons decays with a time constant of ap- 
proximately 10 us or less. Unfortunately, conventional magnet- 
resonance imaging systems are not capable of sampling such ul 
trashort signals. Therefore, the restricted proton pool can be 
depicted only indirectly by exploiting the MT phenomenon _ 

Current methods for determining relaxation parameters of the 
two -pool model including the relative proton density employ both 
the formalism of the coupled Bloch Equations and spectral selec- 
tive RF pulses with the intention to irradiate either the re- 
stricted or the free pool. Following spectral selective radio 





frequency (RF) irradiation, the signal response of the free pool 
is sampled and related to the model parameters . 

In "Magnetic Resonance in Medicine" (Vol. 29, P. 759-766 
(1993)) Henkelman R. M. et al. have derived an equation which re- 
lates the attenuation of the steady state magnetization to the 
resonance offset and power of a continuous wave saturation pulse 
and to some fundamental model parameters including the relative 
size of the restricted pool. The pool parameters can be obtained 
then by fitting this equation to several measurements obtained 
with different RF powers and frequency offsets. So far, this 
method could not be applied in a clinical setting because it ex- 
ceeds the current limit by the specific absorption rate (SAR) and 
because measurement time is impractically long.. Additionally, 
this technique suffers from the constraint of an additional meas- 
urement of the apparent relaxation time Tl. Moreover, prior 
knowledge regarding the lineshape and the transverse relaxation 
time of the restricted pool is required. 

Instead of off -resonant saturation pulses, another method, 
presented by Daniel Gochberg et al. in "Magnetic Resonance in 
Medicine" (Vol. 41, P. 1065-1072 (1999)), employs a train of short 
on-resonant inversion pulses. This approach benefits from a low 
SAR as the inversion pulses are separated by 120 ms. However, it 
suffers from a long measurement time as well. If more than one 
slice is required to be imaged, measurement time scales linearly 
with the number of slices as this method is not capable of true 
multislice imaging. So far, no in vivo results have been gener- 
ated with this method. 

More recent work, presented by John Sled et al. in "Magnetic 
Resonance in Medicine" (Vol . 46, P. 923-931 (2001)) and by Vasily 
Yarnykh in "Magnetic Resonance in Medicine" (Vol . 47, P. 929-939 
(2002)) is also based on the steady state solution of the coupled 
Bloch Equations but regards pulsed RF saturation by incorporating 
numerical calculations. However, these newer methods are still 
not suited for routine use due to long acquisition times and ex- 
tensive calculations . 

It is an object of the present invention to provide a above 
described method which is capable of eliminating the drawbacks of 
the prior art and which enables a fast and low specific absorp- 
tion rate (SAR) method for measuring the macromolecular proton 
density. 
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The present invention relates to magnetic resonance spec- 
troscopy or imaging techniques in which an initial preparation of 
the proton magnetization is performed by two radio frequency 
pulses which causes all, or a portion of, the proton magnetiza- 
tion in the volume of interest to be oriented longitudinally. The 
application of a gradient field between these RF pulses causes 
the longitudinal magnetization to be modulated along the field 
direction. Such a preparation scheme, commonly known as stimu- 
lated echo preparation, does only affect protons of the free mo- 
bile tissue water. The restricted proton pool is not affected, 
because a RF pulse separation time is chosen that is much longer 
than the transverse relaxation time of the restricted proton 
pool. The present invention, therefore, employs above preparation 
scheme such as to label the free water protons and uses these 
protons as an intrinsic indicator to measure the size of the re- 
stricted pool. Following labelling, the concentration and there- 
fore the signal intensity of the labelled protons will decrease 
as the labelled protons "dilute" into the restricted proton pool 
via magnetization transfer. Additionally, the labelled protons 
are also subject to longitudinal relaxation. Probing the labelled 
magnetization by means of a third RF pulse after different mixing 
times yields an indicator decay curve. From this curve the re- 
laxation and dilution effects are separated by a bi -exponential 
analysis and the size of the restricted proton pool is calculated 
according to indicator dilution theory. 

In accordance with a preferred embodiment of the invention, 
an additional RF pulse is applied between the labelling pulses 
and the final readout pulse. This additional pulse is intended to 
change the magnetization in the free proton pool only. Addition- 
ally, a mixing time, i.e. the time between the second and the 
last RF pulse, is used that is twice the time needed for both 
proton pools to restore equilibrium. Then, the size of the re- 
stricted pool is determined from two acquisitions, where the flip 
angle of said RF pulse is 0° in one run and 180° in the other 
run. Said RF pulse can also be a composite pulse that provides a 
constant specific absorption rate independent of the effective 
flip angle. 

Alternatively, in a variation of the method, the flip angle 
of said RF pulse is held constant throughout several runs, where 
the time 1 between said RF and the final readout pulse is varied. 



Such a variation allows to make the mixing time shorter than it 
takes to restore equilibrium. In this fashion, scan time can be 
reduced. 

In an additional variation of the method, in accordance with 
the invention, said RF pulse is applied with a resonance offset 
such as to selectively saturate a part of the restricted pool 
magnetization. Thereby one obtains information about the spatial 
distribution of different spectral components of the restricted 
pool. In accordance with this variation, said off-resonant RF 
pulse can be replaced by a train of off-resonant RF pulses. In 
this manner, a higher degree of saturation and a higher spectral 
selectivity are achieved. 

Further characteristic features of the invention and those 
already mentioned above will be explained in more detail by way 
of the accompanying drawings, wherein Fig. 1 shows a block dia- 
gram of a magnetic resonance imaging system that may be pro- 
grammed to measure the macromolecular proton density; Fig. 2 is a 
simplified model of relaxation in heterogeneous tissues; Fig. 3 
illustrates the principal of the indicator dilution technique to 
measure the volume of distribution or the fractional volume; Fig. 
4 shows a preferred MR pulse sequence for* the determination of 
the macromolecular proton fraction; Fig. 5a-5c show the mapping 
of the macromolecular content in brain tissue with the preferred 
embodiment of the invention; and Fig. 6 illustrates a variation 
of the preferred pulse sequence. 

Below, a technique for measuring the macromolecular *H den- 
sity in relative or absolute terms in a two-step or multi-step 
scan will be described with references to Figs. 1-6. The descrip- 
tion given herein with respect to those figures is for explana- 
tory purpose only and is not intended in any way to limit the 
scope of the invention. 

A system for acquiring the data and generating the images is 
shown in Fig. 1. This system can be a 1.5T or 3T whole-body sys- 
tem from Philips Medical System, Best, The Netherlands, or any 
other suitably equipped MRI system that may be programmed to 
measure the macromolecular proton concentration. As illustrated, 
a magnet 10 generates a static, fundamental magnetic field along 
a z axis 12, in which an object or the body of a patient 14 to be 
examined is situated. The system additionally comprises gradient 
amplifiers 16, gradient coils 18 , transmitter 20, RF power am- 
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plifier 22, RF coils 24 for generating pulse sequences for appli- 
cation to selected slices of the patient's 14 anatomy or the sam- 
ple. The control of the pulse sequence is done by a sequence 
control unit 26, which can be programmed through a scan control 
interface 28. Since software techniques for generating pulse se- 
quences with the characteristics defined below are believed to be 
well-known to those skilled in the art, such pulse generating 
techniques will not be described in further details herein. The 
signal generated by the pulse sequence is received by signal re- 
ceiver 30 and digitized at digitizer 32 for application to an 
arithmetic unit 34 for processing in accordance with the tech- 
nique of the invention. The processed signal is then displayed on 
a display unit 36. Data storage 38 and filming. with a camera de- 
vice 40 may be provided additionally. To synchronize the pulse 
sequence with physiological signals from patient 14, a cardiac 
synchronization unit 42 and a controlling device 44 for respira- 
tory motion may also be provided. 

Now, suppose water protons in tissue to be in two different 
states in regard to their molecular mobility. The compartmentali- 
zation in terms of molecular mobility is illustrated in Fig. 2. 
This figure shows a simplified model of relaxation in heterogene- 
ous tissues. This is commonly know as the two-pool model. The 
pool A corresponds to *H spins in "free" mobile bulk water, 
whereas the second pool B, usually much smaller in size, corre- 
sponds to motional restricted X K spins bound to macromolecules . 
Pool B will be referred to as restricted or macromolecular pool. 
Each pool is characterized by its intrinsic relaxation rates Ri 
and R 2 and by its size Ma° and Mb 0 . Most importantly, there is an 
intermediate to fast exchange of magnetization between both pools 
as indicated by the first order forward and backwards transfer 
rates k f and kb. The fundamental parameter which can be measured 
with the present invention is the size of the macromolecular 
pool. For simplicity the molar fraction f will be used in further 
calculations, which is the relative size of the macromolecular 
pool and which can be written as 
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Since the condition for microscopic reversibility has to be ful- 
filled, Ma°'k £ = M b 0< kb and Equation [1] can also be written as 



As known by those skilled in the art of biomedical engineering, a 
common concept to determine the unknown distribution of volumes 
or fractional sises is to measure the dilution of an indicator. 
The indicator can be of any type, but must be inert and directly 
accessible to a measurement such as in pool A in Fig. 3(a) for 
instance. After application of the indicator into pool A, indi- 
cated by the filled circles/ an exchange process with pool B will 
start. As soon as the dilution of the indicator, caused by diffu- 
sion or other exchange processes, has reached a steady state, 
i.e. the concentration of the indicator is the same in both pools 
(Fig. 3(b)), the relative volume fraction or pool size can be 
calculated from the change in indicator concentration because 



where [c 0 ] is the initial concentration in pool A and [c 8S ] is the 
steady state concentration. 

The present invention uses labelled spin, magnetization as an 
inherent and exchangeable indicator. Spin labelling is achieved 
with a stimulated echo preparation scheme, consisting of two suc- 
cessive RF pulses, preferably with a flip angle of 90° , and a 
gradient field between. For those skilled in the art it is well 
known, that the gradient field causes a phase modulation of the 
transverse magnetization along the field direction. After apply- 
ing the second RF pulse, the phase modulation turns into an modu- 
lation of the longitudinal magnetization. The present invention 
uses a RF pulse separation which is much longer than the trans- 
verse relaxation time of the restricted proton pool. This en- 
sures, that only spins in the free proton pool are labeled even 
if the RF pulses partly saturate the restricted pool . Following 
the generation of this "indicator", subsequent decay of the la- 
belled spins with increasing mixing time is caused primarily by , 
two processes. These are the dilution effect and Tl relaxation. 




[2] . 




[3] 



This process can be modelled by regarding relaxation and magneti- 
zation exchange in the two pool system: 

dM (/) 

--^=-M a (t)(R la +k f )+M b (t)k b [4] 

and r 
dM b (t) 



^ =-M b (t)(R Jb +k b )+M a (t)k f [5] 

where the subscript a and b refer to the free and restricted pool 
and M(t). is the magnetization of the labelled l H spins . 

Solving Equations [4] and [5] with the condition of Mb(t=0) 
= 0 gives a biexponential function for the decay of the labelled 
spins. Only the solution for the free pool is considered here be- 
cause the bound pool will not contribute to the measured signal 
due to its extremely high R 2 : 

M a (t)=M 0 (C 1 expf-A, t)+C 2 exp(-A 2 1)) [ 6 ] 

with 



A 1<2 = 0 . 5 (k b +k f +R la +R lb )±.0 . 5 V(k b +k f +R la +R lb ) 2 -4(R la R lb +R la k b +R lb k f ) 



and 



A y -k h -R lh 
A,-A 2 



A t -A 2 



where Mo is the maximum magnetization available immediately after 
labelling. The rate A* is a "fast" rate which is responsible for 
the quick approach to a steady state between both pools. In con- 
trast, A 2 is a "slow" rate and is roughly comparable to Ri ob- 
tained from a conventional Ti measurement. It is notable that A* 
and A 2 are identical to the rates given by the general solution 
of the coupled Bloch equations after disturbance of the equilib- 
rium (Journal of Magnetic Resonance, Vol. 31, P. 207-229 (1978)). 
However, the constants. Ci and C 2 found here differ from this solu- 
tion. 

In most tissues the backward transfer rate k b is much higher 
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than the longitudinal relaxation rates. At 1.5 T h is approxi- 
mately 16 times higher in white matter, 35 times higher in gray 
matter and 70 times higher in muscle tham any of the longitudinal 
relaxation rates (Magnetic Resonance in Medicine, Vol. 33, P. 
476-482 (1995), Magnetic Resonance in Medicine 35:277 (1996)). 
Therefore, the condition kb»Ria,Rib allows to reduce the constants 
Ci and C 2 to: 




[10] 



and 




til] 



Rewriting Eq. [6] for the measured signal intensity obtained from 
a stimulated echo will result in : 

5(0= ^7^ (/eXp( " A ^ )+eXp( ^0) [12] 

where So is the maximal possible signal intensity that would be 
expected for a mixing time of zero. Thus, f can be calculated 
from a bi -exponential fit to a stimulated echo data set obtained 
with several different mixing times . It is important to note that 
f then is independent of the initial state of both pools. Inef- 
fective, labelling of the free pool will result in a lower dynamic 
range of the decay curve, but this will be reflected only in a 
decrease of So and in the goodness of the fit. The result is also 
insensitive to any incidental irradiation of the restricted pool 
as long as this does not occur during the mixing period. For a 

mixing period much longer than 1/Xi (at 1.5 T this corresponds to 
150 - 200 ms in brain white matter), Equation [12] can be rewrit- 
ten as 

S(0=V^exp(-A a *) [13]. 

When compared to Equation [3] it is obvious that labelled spins 
can be treated as an indicator provided longitudinal relaxation 
is considered. 

Now, consider a two-point pulse sequence such as in FIG. 4, 
which makes use of the dilution effect of phase modulated spins . 
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Fig. 4 shows a preferred MR pulse sequence for the determination 
of the macromolecular proton fraction. This sequence exemplarily 
employs gradients for slice selection (S) , phase encoding (P) , 
frequency encoding (R) , and labelling (M) . The sequence is per- 
formed twice without and with a 180° RF pulse placed in the cen- 
ter of the mixing period. The first run serves as reference scan 
to asses the magnitude of the relaxation term exp(-A, 2 TM) which is 
equal for both runs. The intention of the second run is to intro- 
duce an imbalance in the system. This is achieved with the 180° 
inversion pulse which changes the sign of the phase of the spins 
in the free pool only. This will be followed by a dilution of 
spins with opposite phase to the restricted pool and by a dilu- 
tion of originally labelled spins from the restricted pool to the 
free pool. Immediately before the 180° RF pulse the signal from 
the labelled spins is 

s $ =s 0 j^vL V {-\ 2 m l ) [i4] . 

The signal intensity achieved in the first run is 

S=S'exp{-\ 2 TM 2 ) [15] 

and for the second run with the 180° RF pulse the signal inten- 
sity is 

S 2 --S'-^exp(-A 2 rM 2 )+^ [16] 

where the second term in Equation [16] gives the dilution from 
originally labelled spins in the restricted pool to the free 
pool. In above equations the inversion pulse is assumed to be 
perfect. However, the inversion pulse may also affect the spins 
in the restricted pool which is considered by the parameter p 
that gives the relative saturation of the restricted pool; total 
saturation will result in a p value of 0 while no saturation will 
be indicated by a value of 1. Usually, the restricted pool is ex- 
pected to be smaller then the free pool (f<l) . Then, if the sig- 
nal intensity is obtained from a magnitude . image the molar 
fraction is given by 

/= <^7iJ 1171 • 



If R2b is known or can be estimated, p can be obtained from nu- 
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merical simulations of the effect of the inversion pulse on the 
bound pool. Equation [17] can be implemented to calculate pixel- 
by-pixel parameter images showing the molar fraction. To measure 
the absolute macromolecular proton concentration, the proton con- 
centration of the free water pool has to be determined. This can 
be achieved, for instance, by estimating the size of the long T 2 - 
component from an additional multi-echo experiment and by scaling 
it to a reference water sample of known temperature . 

The pulse sequence of Fig. 4 was implemented on a Intera™ 
1.5T whole-body scanner. 

The pulse sequence and method was validated with phantoms 
containing known concentrations of agar gel and bovine serum al- 
bumin. Additionally, the pulse sequence was used to evaluate the 
macromolecular content in brain tissue of several volunteers. 
Fig. 5 illustrates an application of the method of mapping the 
relative proton density of the macromolecular pool in a patient 
suffering from multiple sclerosis (MS). Figs. 5(a) and 5(b) are 
the images corresponding to the two acquisitions performed with 
and without the inversion pulse a*, respectively. Fig. 5(c) shows 
the macromolecular fraction computed pixel-by-pixel according to 
Equation [17] . In brain tissue the macromolecular pool is associ- 
ated to the myelin lipids and proteins and is therefore expected 
to scale with myelin density. It can be nicely appreciated that 
white matter exhibits a higher macromolecular proton density than 
the gray matter does. MS plaques show a reduction in the macro- 
molecular proton density. Of course, the macromolecular proton 
density in other tissues such as myocardial tissue or cartilage 
can be similarly determined using the techniques of the inven- 
tion. ^ 

Although exemplary embodiments of the invention have been 
described in detail above, those skilled in the art will appreci- 
ate that many additional modifications are possible in the exem- 
plary embodiment without materially departing from the novel 
teachings and advantages of the invention. The pulse sequence il- 
lustrated in Fig. 4 represents only some of many possibilities 
for measuring the macromolecular proton density from the dilution 
of phase modulated spins. In this sequence, a RF pulse (Ob) is 
applied to disturb the equilibrium of the labelled spins by ma- 
nipulating the spins of the free proton pool. Alternatively dis- 
turbance " of the equilibrium can be achieved with a sequences like 
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in Fig. 4, but which uses single or trains with flip angles 
a 3 ,i....<X3,n of off -resonant RF pulses instead of the third RF 
pulse of flip angle oc 3 in order to saturate the spins in the re- 
stricted proton pool. Such variation of the preferred pulse se- 
quence is illustrated in Fig. 6. Also, this principle can 
likewise be applied to sequences like in Fig. 4, but where the 
flip angle of the inversion pulse deviates from 180° , or to a 
scheme were Ta or X* are varied over successive scans. Accord- 
ingly, all such modifications are intended to be included within 
the scope of this invention as defined in the following claims. 




Claims : 



1. A method of measuring the dilution of phase modulated spins, 
from which the macromolecular proton concentration involved in 
magnetization transfer can be calculated, from two or more scans 
of an object with a magnetic resonance imaging device, comprising 
the steps of: 

- applying a first RF pulse of flip angle Ofc at a first time so 
as to generate a transverse magnetization in said object; 

- applying a first magnetic field gradient along a predetermined 
direction in said object so as to produce a phase modulation of 
*H spins along the direction of the gradient; 

- applying a second RF pulse of flip angle a 2 at a second time Ti 
seconds after said first time so as to flip the transverse mag- 
netization into the longitudinal plane; 

- applying a third RF pulse of flip angle Oa at a third time T1+T2 
seconds after said first time so as to manipulate the longitudi- 
nal stored magnetization or as to leave the longitudinal stored 
magnetization unaffected; 

- applying a /fourth RF pulse of flip angle cu at a fourth time 
Xi 4-t 2 +t 3 seconds after said first time so as to flip the longitudi- 
nal stored magnetization into the transverse plane; 

- applying a second magnetic field gradient along the same prede- 
termined direction as said first magnetic field gradient; 

- detecting a stimulated echo at a fifth time. 

2. A method according to claim 1, wherein the spins of 13 C, 14 N, 
19 F, 23 Na, 31 P, or 35 C1 nuclei are phase modulated. 

3. A method according to claim 1 or 2, wherein said third RF 
pulse is a composite RF pulse. 

4. A method according to claims 1 to 3, wherein said third RF 
pulse is applied with a resonance offset so as to saturate partly 
or fully the magnetization associated to the macromolecular pool. 

5. A method according to any one of claims 1 to 4, wherein said 
third time is altered over several subsequent scans. 

6. A method according to any one of claims 1 to 5, wherein said 
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pulse sequence is a pulse sequence for multislice imaging. 

7 A method according to any one of claims 1 to 6, wherein said 
Pulse sequence incorporates an additional encoding gradient for 
3D imaging. 

8. A method according to any one of claims 1 to 7, wherein said 
fourth RF pulse is replaced by a train of rf pulses with a flip 
angle < 90° in order to acquire more than one line in k-space per 
repetition. *F*ee per 

9 a method according to any one of claims 1 to 7, wherein said 
stimulated echo is sampled with a multi-shot or single-shot echo 
planar imaging technique. 

10. A method according to any one of claims 1 to 9, wherein the 
object is a patient, and wherein the first RF pulse is synchro- 
nized using electrocardiographic gating or peripheral pulse gat- 

11. A method according to any one of claims 1 to 10, wherein the 
object is a patient, comprising controlling the respiratory mo- 
tion of the patient during application of the pulse sequence. 

12 A method according to any one of claims 1 to 11, comprising 
determining a longitudinal relaxation rate from the sampled data, 
whereas a spin echo is sampled additionally at a time 2* after 
said first time. 

13. A method according to any of claims 1 to 12, wherein said ob- 
ject is the brain tissue of a patient and the macromolecular pro- 
ton concentration represents myelin density of said patient. 

14 A method according to any of claims 1 to 13, wherein said ob- 
ject is the myocardium of a patient and the macromolecular proton 
concentration reflects fiber density and structure and therefore 
•tissue quality. 

15. A method according to any of claims 1 to 14, wherein the mac- 
romolecular proton pool is a contrast agent administered to the 
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TritlLTf T maCr ° raoleCular P~ton density reflects the concen- 
tration of the contrast agent. 

16 A magnetic resonance imaging device comprising a magnet (10, 

coils (18, which apply gradient pulses to said object (14) rf 
coils ,24, which apply RF pulses to said object <L,, driving 

™7u ( 22, . WhiCh COils <") -d RK 

coils (24,, receding circuitry (30, which receives a signal from 

IT1T C \ 114) in Said ma9netiC fi6ld UP ° n a ^i-tion^f said 
gradient pulses and rf pulses, an arithmetic. unit (34), a displav 
device (36, for displaying said received and processed signals 
and a sequence control device (26, which Controls said rf coils 
(24) to generate and apply a first RF pulse of flip angle a, at a 

object M4Y 0 r t0 gSnerate a averse magnetization in said 
object (14, and a second RF pulse of flip angle cc 2 at a second 
txme Tl seconds after said first time so as to flip the trans- 
verse magnetization into the longitudinal plane and a third RF 
Pulse of flip angle oc 3 at a third time t 1+ t 2 seconds after said 
first time and a fourth rf pulse of flip angle a. at a fourth 
time t 1+ t 2+ t 3 seconds after said first time so as to flip the Ion 
gitudinal stored magnetization into the transverse plane and 
which sequence control device (26, controls said gradient coils 
(18, to generate and apply first and second magnetic field gradi- 
ents along a predetermined direction in said object (14, and 
which generates an image of a stimulated echo detected by said 
receiving circuitry (30, at a fifth time. 

17 An imaging device according to claim 16, wherein said arith- 

18. An imaging device according to claim 16 or 17, wherein said 

RF eq puL C e e ZT* d6ViCe (26> ^ ^0 perform said third 

RF pulse with a resonance offset so as to saturate partly or 
fully the magnetization associated to the macromolecular pool. 

19. An imaging device according to any one of claims 16 to 18, 
wherein said sequence control device (26, is programmed for al- 
tering said third time over several subsequent scans. 
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20. An imaging device according to any one of claims 16 to 19, 
wherein said sequence control device (26) controls the gradient 
amplifiers (16) such as to sample a signal with a digitizer (32) 
according to the multi-shot or single-shot echo planar imaging 
technique. 

21. An imaging device according to any one of claims 16 to 20, 
characterized in a synchronization unit (42) connected with a de- 
vice for measuring the electrocardiographic activity of a patient 
for synchronisation of the first RF pulse with the electrocardio- 
graphic activity of the patient. 

22. An imaging device according to any one of claims 16 to 21, 
characterized in a controlling device (44) for controlling a res- 
piratory motion of a patient during application of the pulse se- 
quence . 

23. An imaging device according to any one of claims 16 to 22, 
wherein said sequence control device (26) is additionally pro- 
grammed to sample a spin echo at a time 2Ti after said first 
time, in order to determine the longitudinal relaxation rate. 
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ABSTRACT 

A method and system for measuring the macromolecular proton con- 
centration involved in magnetization transfer comprising the 
steps of 

- applying a first RF pulse of flip angle Ok at a first time, 

- applying a first magnetic field gradient along a predetermined 
direction, 

- applying a second RF pulse of flip angle a 2 at a second time Ti 
seconds after said first time, 

- applying a third RF pulse of flip angle Ota at a third time Ti+t 2 
seconds after said first time , 

- applying a fourth RF pulse of flip angle ou a,t a fourth time 
T1+T2+T3 seconds after said first time, 

- applying a second magnetic field gradient along the same prede- 
termined direction as said first magnetic field gradient, and 

- detecting a stimulated echo at a fifth time. 
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